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ABSTRACT 

Large surveys for Lya emitting (LAE) galaxies have been proposed as a new 
method for measuring clustering of the galaxy population at high redshift with 
the goal of determining cosmological parameters. However, Lya radiative trans- 
fer effects may modify the observed clustering of LAE galaxies in a way that 
mimics gravitational effects, potentially reducing the precision of cosmological 
constraints. For example, the effect of the linear redshift-space distortion on 
the power spectrum of LAE galaxies is potentially degenerate with Lya radia- 
tive transfer effects owing to the dependence of observed flux on intergalactic 
medium velocity gradients. In this paper, we show that the three-point function 
(bispectrum) can distinguish between gravitational and non-gravitational effects, 
and thus breaks these degeneracies, making it possible to recover cosmological 
parameters from LAE galaxy surveys. Constraints on the angular diameter dis- 
tance and the Hubble expansion rate can also be improved by combining power 
spectrum and bispectrum measurements. 

Key words: galaxies: high-redshift - cosmology: theory - large-scale structure 
of universe 



1 INTRODUCTION 

For the past three decades galaxy redshift surveys have 
served as the traditional method for constraining cosmo- 
logical parameters such as the matter density of the uni- 
verse and the equation of state of dark energy, by measur- 
ing the clustering of galaxies. These have been restricted 
to 2 < 1 due to the increasingly fainter galaxy magni- 
tudes and larger required cosmic volumes, which render 
spectroscopy of large numbers of photometrically selected 
early type galaxies plausible only at such low redshifts. 
Recently the WiggleZ collaboration has pushed galaxy 
clustering work to z ~ 1 using emission lines from star- 
forming galaxies (Bla ke et al.|2011a|b|c 20121. 

Lya emitting (LAE) galaxies are detectable out to 



* E-mail: bgreig@student.unimelb.edu.au (BG) 
f komatsu@mpa-garching.mpg.de (EK) 
I swyithe@unimelb.edu.au (JSBW) 



high redshift ( 


lye et al. 2006 Kashikawa et al. 2006 Lehn- 


|ert et al.|2010 


|Ouchi et al.|2010|), due to their strong line 



emission. Indeed over the previous few years, the number 
of detected LAE sources has steadily grown and the sam- 
ple sizes of LAE galaxies have reached sufficient size for 



clustering studies ( Gawiser et al.|2007||Kovac et al.|2007| 
Orsi et af |2008] |Guaita et al.|2010||Ouchi et al.||2010[ ). 



While the existing samples of LAE galaxies are still 
too small for cosmological purposes, the rate of detection 
of these LAE galaxies will significantly improve with the 
upcoming Hobby-Eberly Telescope Dark Energy Experi- 
ment (HETDEX, |Hill et al.|2004l [2008] ) , whose aim is to 
spectroscopically measure the redshifts of 800 000 LAE 
galaxies in the redshift range 1.9 < z < 3.5 (Hill et al 



2004 2008 1 with the total sky coverage of 420 square de- 
grees and the total volume coverage of 10 Gpc 3 . This sur- 
vey is specifically designed to use the clustering of LAE 
galaxies to make the precise measurement of the distance 



© 0000 R AS 



2 B. Greig et al. 



scales, both the angular diameter distance {Da) and the 
Hubble rate (H), as a function of z out to z ~ 3. 

In order for us to use the clustering of LAE galaxies 
to measure the distance scales, we must understand how 
the clustering of LAE galaxies is related to the underlying 



matter distribution. Simulations by Zheng et al. (2010 



2011 1 and Laursen et al. (2011 \ have investigated the ra- 



diative transfer effects on the Lya emission of LAE galax- 
ies both within the circumgalactic environment around 
the halo and from the resonant scattering of diffuse neu- 
tral hydrogen in the intergalactic medium (IGM). Of 
particular interest is the clustering of the LAE galax- 
ies: Zheng et al. (20111 find that line-of-sight gradients 



in the peculiar velocity of LAE galaxies could lead to an 
observed reduction in the line-of-sight clustering ampli- 
tude of the galaxies, counteracting the strength of the 
typical Kaiser effect (Kai ser|1987[ ) caused by gravitation. 
Conversely, the clustering of LAE galaxies transverse to 
the line-of-sight is found to be significantly boosted by 
the Lya radiative transfer effects. In addition to the lo- 
cal effects of peculiar velocity gradients, the transmission 
of the Lya emission line of LAE galaxies through the dif- 
fuse IGM could also be affected by fluctuations in the 
UV ionizing background and to changes in the neutral 
hydrogen fraction associated with changes in the density 
around the local environment. 



To understand these effects, Zheng et al. (2011 1 and 



Wyithe & Dijkstra (20111 have derived analytic mod- 



els to describe the observed modifications of the power 
spectrum of LAE galaxies. Both derive quantities that 
directly relate to the non-gravitational effects expected 
from the Lya radiative effects. Wyithe & Dijkstra (2011 1 



use this model to study the expected recovery of both cos- 
mological and Lya radiative transfer parameters from a 
survey corresponding to HETDEX. They find that some 
cosmological parameters derived only from the power 
spectrum are degenerate with the Lya radiative transfer 
effects, and that this has direct consequences for the ac- 
curacy with which cosmological parameters can be recov- 
ered from the LAE galaxy power spectrum. Prior knowl- 
edge of the magnitude of the radiative transfer effects can 
improve the recovery of the cosmological constraints. 

In this paper we further investigate the effects of 
non-gravitational LAE clustering on the recovery of cos- 
mological parameters. We extend and improve the linear 



theory work of Wyithe & Dijkstra (20111 by including 



the three-point correlation function (bispectrum) and by 
combining with the power spectrum, to break the first or- 
der degeneracies of the Lya radiative transfer effects and 
cosmological parameters. To calculate the bispectrum, we 
use a next-to-leading order Eulerian perturbation theory 
approach ( Bernardeau et al.|2002 and references within) , 
and derive expressions valid into the mildly non-linear 
regime. We also derive higher-order expressions for the 
Lya radiative transfer effects, including the higher-order 
effects of redshift-space distortions. We then study how 
well a joint analysis of the power spectrum and the bis- 
pectrum can break the cosmological and radiative trans- 
fer degeneracies. We provide the expected constraints 
on cosmological parameters through the application of 
Fisher matrices, with specific reference to the HETDEX 
survey. 



This paper is set out as follows. In Section [2] we out- 
line the degeneracies between cosmological and Lya ra- 
diative transfer parameters, and in Section [3] perform a 
Fisher matrix analysis of the LAE galaxy power spec- 
trum. In Section [4] we outline and describe existing Eu- 
lerian perturbation theory expressions and provide the 
derivation of higher-order corrections for the Lya radia- 
tive transfer effects, in order to construct both a bispec- 
trum and a reduced bispectrum model. In Sections [5] and 
[6] we perform Fisher matrix analyses of the reduced bis- 
pectrum alone, and a combined power spectrum and bis- 
pectrum in order to provide cosmological parameter es- 
timates. We finish with our summary and final remarks 
in Section In our numerical calculations we consider 
the standard set of cosmological parameters ( |Komats"u" 
etllpHlj ), with D. m = 0.275, ft A = 0.725, Q b = 0.0458, 
n s = 0.968, h = 0.702 and a 8 = 0.816. 



2 CLUSTERING OF Lya EMITTERS 

In this section we summarize the linear theory clustering 
of LAE galaxies. For a galaxy redshift survey, one can 
write a simple expression relating the power spectrum 
of galaxies to the underlying matter distribution. How- 



ever, for LAE galaxies, Zheng et al. ( 2011 \ and Wyithe & 
Dijkstra (20111 show that Lya radiative transfer effects 



modify this relationship. 



2.1 Galaxy power spectrum 

In linear theory, the power spectrum of galaxies in red- 
shift space is given by 



(1) 



where b\ is the linear galaxy bias, since galaxies are biased 



tracers of the matter density field (Kaiser 1984), / = 
din D(a) / din a is the growth rate of structure, D(a) is the 
linear growth factor, Pi (ft) is the power spectrum of the 
linear density fluctuations, and fj, is the cosine of the angle 
between the line-of-sight vector z and the wavevector k, 
i.e., fi = k ■ z/k. 

Of particular interest for cosmology is the recov- 
ery of the growth rate of structure, /, which can be 
parametrized as f — fi^. The value that 7 takes can dis- 
tinguish between cosmological models described by gen- 



eral relativity or other gravitational descriptions (Lin- 
der 2005). However, the galaxy bias, 61, and the growth 
rate of structure, /, are degenerate in the above model. 
Hence, what is actually measurable from a galaxy redshift 
survey is the linear redshift-space distortion parameter, 
ft — With this parameter, the above expression be- 



P gal (fc) = &l(l + /3M 2 ) 2 Pi(A0- 



2.2 LAE power spectrum 



(2) 



Lya radiative transfer effects potentially modify the 
power spectrum of LAE galaxies away from the standard 
galaxy description (Equation [TJ . This arises because the 
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observed number density of Lya galaxies at a fixed ob- 
served flux depends on how many Lya photons escape 
to observers. Thus the observed density of galaxies can 
be modified by the local environment nearby to the Lya 
galaxies. 

The Lya optical depth through the IGM depends on 
the local density (p), the ionizing background strength 
(r), and also on the strength of the local peculiar velocity 
gradient along the line-of-sight (dv z /dr), 



pyo.7 di> z 

dr 



-2-0.7(7-1) 



dr 



(3) 



where 7 is the polytropic index, used to relate the tem- 
perature to the underlying density field as T oc p 7_1 with 
7 = 1.4 (|Hui fc Gnedhi|l997l. 



In this section we shall describe the effect of Lya 
transmission fluctuations on the observed clustering sig- 



nal of LAEs, following |Wyithe fc Dijkstra| ( |201l| ). Modifi- 
cations to the intrinsic Lya luminosity by the IGM induce 
a change in number counts of LAE galaxies. The number 
density of LAE sources, nL yo , that are observed above 
some observational flux threshold, Fo, can be expressed 
relative to the average, ri,Lya(> Lo, po, To, S(x)), as 

n Ly a(> F ) = n Lya (> L ,p ,r ,8(x)) x [I + 5 g (x)], (4) 

where S(x) is the large-scale matter-density perturbation, 
and 8 g (x) is the perturbation in the galaxy number den- 
sity. Now the average number density of LAE galaxies 
depends on the fluctuations due to the non-gravitational 
Lya radiative transfer effects; namely, the local density 
in the LAE environment, p, ionizing background, V, and 
peculiar velocity gradients, di; z /dr. 

We first obtain the expression for the mean num- 
ber density of observed LAE galaxies. Taylor expanding 
about the three radiative transfer effects, we obtain 



n Lya (> L ,p ,T ,S(x)) = n$ a ( 



1 + n 



(i) 

Lya 



(5) 



where n^ a is just the mean number of LAE galaxies, 

and n^ya ' s the first-order Taylor-expanded expression 
evaluated around their mean quantity, 



l (0) 
l Lya 



-(r-r 



<9nLyo 



dv z 
d(ar cc 



dT 



H 



+ (P 
F ,r 



Po) 



On 



Lya 



dp 



d 



dv z 



d(ar com ) 



F(),P0 



(6) 



Fo.PO 



Here, H is the Hubble rate, and the line-of-sight velocity 
gradient is taken with respect to the comoving distance, 



We rewrite this expression for n\^! a as 



Lya 



Lya 



= $r(x)Cr + S p (x)C p + S v (x)C v . 



(7) 



The constants Cr, C p , and C„, are defined in Appendix 
[A} and capture the distinct physical effects caused by 
changes to the local environment corresponding to either 
changes to the ionizing background, the density, or the 
velocity gradient along the line-of-sight, respectively. The 
other quantities are defined as Sr = (T — ro)/Lo, S p = 
(p - pa) /pa, and S v = (Ha)~ 1 dv z /dr coln . 

Next, the number density of LAE galaxies observed 



above a flux limit can be written as, 

n L ya(> -Fo) 



(8) 



where we have substituted Equation [5] into Equation [4] 
Rewriting this expression as fluctuations in the number 
of LAE galaxies relative to the mean number of galaxies 
expected without Lya radiative transfer effects, nQ a , we 
obtain 



Shya{x) 



riLya(> F ) 

-(0) 
n Lya 



\l+d g (x)] [l + 4^-1.(9) 



Expanding Equation [9] and taking the Fourier transform, 
one finds, to the first order, 

5 L ya, s (fc) = [61 (1 + Cr) + C P + fp 2 (l - C v )] S(k), (10) 

where we have used the linear-theory predictions: 
S r (k) = biC r S(k), S p {k) = S(k), and S v {k) = 
—fp 2 C v S(k). Here, we have implicitly performed the lin- 
ear redshift-space transformation (see Appendix [5]) . 

In this work we define the redshift-space LAE galaxy 
power spectrum as 

<<5Lya, S (fel)<5Lya, S (fc 2 )) = (2tt) 3 P L ya, 3 {k) 5° (fcl + fc 2 ) , (II) 

where PLya,s(fc) is given by 

P L ya,s(fc) = [6l (1 + Cr) + C p + //i 2 (l - C„)] 2 P L (k), (12) 

and Pl (k) is the linear real-space matter power spectrum. 
Relative to Equation [l] the addition of Lya radiative 
transfer effects can lead to changes in the amplitude of 
the measured power spectrum^ 

Equation |12| cont ains the main contributing terms 



of Zheng et al. (20111. However, we do not include the 



transverse line-of-sight velocity-gradient or the density- 



gradient (as provided by Zheng et al. 20111. As shown 
in Zheng et al. (20111, the effect of the density-gradient 



adds additional scale-dependant terms to the expression 
for the clustering of LAE galaxies on small scales. In 
this work, we are working at much larger scales, and so 
can ignore this scale-dependence and allow the density- 
gradient terms to be absorbed into the existing parame- 
ters of Equation 1 12| 

The inclusion of Lya radiative transfer effects in- 
troduces degeneracies between the cosmological parame- 
ters and Lya radiative transfer parameters. In particular, 
from Equation |12| we note the degeneracy between the 



1 In Equation ^2] we do not consider the scale dependence of 
the ionizing background fluctuations, which is the major dif- 
ference between our expression and the expression in Wyithc 
|fc Dijkstra| ( |20lH . The ionizing background fluctuations are 
expected to be scale dependant, important on large scales (set 
by the mean free path of the ionizing photons) and becoming 
negligible on small scales. However, including the scale depen- 
dence associated with the ionizing background fluctuations in- 
creases the model complexity, providing additional model de- 
generacies. We feel that this simplification is justified since 
the transmission models investigated by |Wyithe fc Dijkstra| 
poTT} find the magnitude of ionizing background fluctuations 
(Cr) to be small compared to the other two radiative transfer 
effects. Hence, while ignoring the scale dependence is a simpli- 
fication, the overall impact of removing this scale dependence 
should be minor. 
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growth rate of structure, /, and the line-of-sight pecu- 
liar velocity radiative transfer effect, C v . Additionally, 
the galaxy bias, bi, is degenerate with the local environ- 
ment density, C p , and the fluctuations in the ionizing 
background, Cr. 

To simplify the expression, consider the following re- 
definition of Equation |12| 

fW(fc) = b( [l + /3 M 2 (l-a)] 2 Pi(fc), (13) 
where 

bj^bj+Cp + bjCr- (14) 

This includes the large-scale effects of density and ioniz- 
ing background which modify the observed galaxy clus- 
tering, and P which is a modified linear redshift-space dis- 
tortion parameter corresponding to the modified galaxy 
bias, p = f/b\. 

Now, the problem is clear: while Equation[l3]has the 
same structure as Equation [2] the meaning of each pa- 
rameter is different. The correspondence is bi — > b\ and 
j3 — > /3(1 — C v ), which shows the parameter degeneracy. 
In the next section we illustrate the resulting effect of ra- 
diative transfer parameters on the potential cosmological 
constraints. 



3 COSMOLOGICAL CONSTRAINTS BASED 
ON THE LINEAR LAE GALAXY POWER 
SPECTRUM 

To generate the expected constraints on cosmological pa- 
rameters for a given survey configuration, wc calculate 
the Fisher matrix which, for the galaxy power spectrum, 
can be written as (e.g., Seo &: Eisenstein|2003 1 



k 2 dk 
2tt 2 



ainP, 



Lya ,s 



<91nP, 



Lya,s 



d6i 



where w(k,n) is the weight given by 
w(fc,/i) = 



n g Pr ya , s {k, At) 

oFhva.s (&, /^) 



V„ 



w(k,fi)dfi, (15) 



(16) 



Here, Qi is the parameter set of our ith dimensional 
model, n g is the number density of LAE sources, and 
Vsurvey is the volume of the redshift survey. 

We focus our attention on a survey like HETDEX, 
for which we assume the linear galaxy bias to be bi — 2.2. 
We generate constraints assuming measurement of one 
redshift bin at the midpoint of the HETDEX redshift 
range, z min = 1.9 and ; 



= 3.5. At this redshift we 



have / = 0.972 for the growth rate of structure. We as- 
sume HETDEX will detect 800 000 LAE galaxies in a 
total survey area of 420 sq. deg. We restrict our analy- 
sis to the weakly non-linear regime, selecting a maximum 
wavenumber, fc max = 0.3/iMpc" 1 . 

To generate the l-cr constraints for our cosmological 
parameters, we construct the one-dimensional maximum 
likelihood distribution from the Fisher matrix assuming 
a Gaussian distribution. The likelihood for the ith model 



Parameter 


Marginalization 


PS 






l-cr (per cent) 


/3 


ln(.4) 


0.0091 (2.06) 





ln{A),ln(D A ),ln(H) 


0.0213 (4.82) 


\n(D A ) 


ln(A),/3,ln(J3") 


0.0110 (1.10) 


HH) 


ln(A),/3,ln(D A ) 


0.0132 (1.32) 



Table 1. The l-cr constraints for the linear redshift-space dis- 
tortion parameter, /3, the angular diameter distance, lnjfl^), 
and the Hubble rate, ln(H), for a galaxy redshift survey with 
HETDEX-like survey parameters. The other model parame- 
ters are marginalized over, but no Lya radiative transfer ef- 
fects are included, i.e., the power spectrum is given by Equa- 
tion rj] 



parameter is 



Jzf (xi) = exp 



1-2 



(17) 



where Xi = (xi — x ) is defined to be the cosmological 
parameter value, Xi, subtracted by its fiducial value, x - 
F is the full Fisher matrix of the n parameter model, 
and F the reduced Fisher matrix of the n — 1 parameter 
model with the ith row and column removed. 

To generate the two-dimensional joint constraints on 
any two parameters, we use 



J£{xi,Xj) — exp ■ 



-2 

i 



F t i 



-2 



Fik(F k iT 1 F li 

k,l^i i 



+S1 F. 



Z F ik (F M ) 



L Fu 



(18) 



which contains the cross term which determines the cor- 
relation between the two parameters being considered. 
Here, F is the full Fisher matrix of the n parameter 
model, and F is the reduced Fisher matrix of the n — 2 
parameter model, with the ith and jth rows and columns 
removed. 

To investigate the degeneracies due to Lya radiative 
transfer parameters, we consider recovery of cosmological 
parameters from three power spectra: 

1. The galaxy power spectrum given by Equation [2] 

2. A fiducial LAE power spectrum given by Equa- 
tion [13] with the fiducial values of radiative transfer pa- 
rameters set to vanish, i.e., Cr — C p = C v = (but these 
radiative transfer parameters are marginalized over) , and 

3. A LAE power spectrum given by Equation |13| with 
the fiducial values of radiative transfer parameters set to 
some indicative values. 



3.1 Galaxy power spectrum 

For the galaxy power spectrum given by Equation [2] the 
galaxy bias is completely degenerate with the amplitude 



© 0000 RAS, MNRAS 000, 000-000 



Cosmology from clustering of Lya galaxies 5 



Lfi 

o 
o 



LO 

q 
d 
J_ 

o 



1-0.1 




-0.05 



-0.05 

ln(M 



0.05 



0.1 



Figure 1. Two-dimensional marginalized joint distribution for the two cosmological distance scales: the angular diameter distance 
(Da) ana the Hubble rate (H). (Left) a typical galaxy redshift survey (no Lya effects; no marginalization over C v ), (middle) 
a fiducial LAE galaxy redshift survey (the fiducial values of the Lya radiative transfer parameters set to vanish; marginalized 
over C v ) and (right) a LAE galaxy redshift survey including the first-order Lya radiative transfer effects given by Cr = 0.05, 
C p = —0.39, and C v = 0.11. The error ellipse is slightly bigger for this case because the effective bias of LAE galaxies, b\ = 1.9, is 
about 15% smaller than the fiducial value, b\ = 2.2, reducing the amplitude of the power spectrum relative to the shot noise. The 
solid and dashed curves show the 1- and 2-cr constraints generated from the likelihood distribution, respectively. Scale selected to 
aid comparison with Figure [7] 



Parameter 


Marginalization 


No Priors on C v 
1-cr (per cent) 


Perfect knowledge 
l-<7 (per cent) 


a Cv = 0.01 
l-tj (per cent) 


CT C„ = 0.1 
1-iT (per cent) 


cr Cv = 0.5 
1-cr (per cent) 


P 
P 

ln(D A ) 
ln(H) 


ln(A),C„ 
ln(A),ln(D A ),\n(H),C v 
ln(A)J,ln(H),C v 
ln(A)J,\n(D A ),C v 


0.0110 (1.10) 
0.0132 (1.32) 


0.0091 (2.06) 
0.0213 (4.82) 
0.0110 (1.10) 
0.0132 (1.32) 


0.0101 (2.29) 
0.0218 (4.93) 
0.0110 (1.10) 
0.0132 (1.32) 


0.0451 (10.21) 
0.0491 (11.11) 
0.0110 (1.10) 
0.0132 (1.32) 


0.2211 (50.04) 
0.2220 (50.24) 
0.0110 (1.10) 
0.0132 (1.32) 



Table 2. The 1-tr constraints for p, ln(D^), and ln(H) for our fiducial LAE galaxy redshift survey (the fiducial values of the 
Lya radiative transfer parameters set to vanish) marginalized over remaining model parameters shown in the second column. We 
compare varying priors added to the radiative transfer parameter, C v , which suffers from a large degeneracy with p. The power 
spectrum model is given by Equation |13| 



of the power spectrum (as). Furthermore, one cannot di- 
rectly measure the growth rate of structure, /, but only 
the parameter /3. We show later that by considering the 
bispectrum, one can directly probe /. 

The cosmological parameters we determine from this 
model are therefore the overall amplitude [ln(A)], the lin- 
ear redshift-space distortion parameter [/}], and the two 
distance measurements given by the angular diameter 
distance [ln(D^)] and the Hubble rate [ln(//)]. The in- 
formation on the galaxy bias is factored into the linear 
rcdshift-space distortion parameter and we redefine the 
amplitude to include galaxy bias. The constrains for j3 
as well as ln(DA) and \n(H) represent the best case sce- 
nario for how accurately we can recover the cosmological 
parameters from a galaxy power spectrum analysis. 

In Table [I] we provide the 1-a constraints on /3, and 
the two distance scales, Da and H. For a galaxy red- 
shift survey with HETDEX-like survey parameters, the 
expected uncertainty on the linear redshift-space distor- 
tion, /3, is 0.021 (4.8 per cent), on the angular diameter 
distance it is 1.1 per cent, and on the Hubble rate it is 
1.3 per cent. Our distance constraints for the above model 



are consistent with the results of Shoji et al. (20091. 

Of particular interest for cosmological analyses is 
the two-dimensional joint constraints on the two distance 
measures, ln(DA) and ln(H). In the left panel of Figure 



[T] we show the 1-cr and 2-a joint constraints on \xi(Da) 
and ln(_ff), which give the baseline for comparison with 
the recovery of the cosmological distance parameters for 
Da and H for the remainder of this work. 



3.2 Fiducial LAE galaxy power spectrum 

We model the LAE galaxy power spectrum using Equa- 
tion |13| In our fiducial case, we set all Lya radiative 
transfer coefficients to zero, i.e., 6i = 6i, j3 = /3, and 
C v — 0. Although we set the Lya radiative transfer ef- 
fects to zero, we still marginalize over the possible exis- 
tence of C v in this case. 

The amplitude of the power spectrum is completely 
degenerate with the modified galaxy bias, and so we can 
redefine the amplitude to include 61. The LAE galaxy 
power spectrum includes C v , the Lya radiative transfer 
effect associated with the line-of-sight peculiar velocity 
gradient. Hence, for the LAE galaxy power spectrum, the 
model contains 5 parameters, ln(A), /3, C v , Iti(Da), and 
\n(H). The linear redshift-space distortion parameter, f3, 
and the radiative transfer effect, C v , are completely de- 
generate; however, with the addition of priors on Cv , one 
can break the degeneracy and impr ove the constraints 
on the linear distortion parameter, $ ( Wyithe & Dijkstra 



2011 1 
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Figure 2. One-dimensional marginalized likelihood distributions for the linear redshift-space distortion parameter (/3) generated 
from; (left) the fiducial LAE galaxy power spectrum (no Lyce effects, but marginalized over C„), and (right) a LAE galaxy power 
spectrum including Lyo radiative transfer effects (Cr = 0.05, C p = —0.39, and C v = O.lf). The resultant offset in ft in the 
LAE galaxy power spectrum (right panel) is due to the modified bias because of the included Lya radiative transfer effects. The 
various curves denote different priors added to C v and are as follows, black solid: perfect knowledge (txc„ = 0.0001), grey dotted: 
oq v = 0.01, grey dot-dashed: cq = 0.1, grey dashed: ac = 0.5, and grey solid: no priors added. In the right panel, the black 
dashed offset curve is the comparison to the case corresponding to perfect knowledge on C v from the fiducial LAE galaxy power 
spectrum (black solid curve, left panel). 



Parameter 


Marginalization 


No Priors on C v 
l-o~ (per cent) 


Perfect knowledge 
1-cr (per cent) 


a Cv = 0.01 
1-cr (per cent) 


<?C V = 0-1 
1-cr (per cent) 


(T C „ = 0.5 
1-cr (per cent) 




MD A ) 
\n(H) 


\n(A),\n(D A ),ln(H),C v 
ln(A)J,ln(H),C v 
ln(A)J,ln(D A ),C v 


0.0128 (1.28) 
0.0151 (1.51) 


0.0121 (2.39) 
0.0278 (5.49) 
0.0128 (1.28) 
0.0151 (1.51) 


0.0134 (2.65) 
0.0283 (5.59) 
0.0128 (1.28) 
0.0151 (1.51) 


0.0582 (11.49) 
0.0633 (12.50) 
0.0128 (1.28) 
0.0151 (1.51) 


0.2847 (56.23) 
0.2858 (56.45) 
0.0128 (1.28) 
0.0151 (1.51) 



Table 3. Same as Table [T] but for the LAE galaxy redshift survey including first order Lyct radiative transfer effects, Cr = 0.05, 
C p = -0.39 and C„ = 0.11. 



In Table [2] we provide the resulting 1-cr constraints 
on the linear redshift-space distortion parameter, /3, as 
well as the distance constraints, marginalized over the 
remaining model parameters including C v . The columns 
from left to right in Table [2] consider priors added to 
C v ; essentially perfect knowledge of C v , ac v = O.OOOlQ 
o~c v = 0.01, crc v = 0.1, and ac v = 0.5. The inclusion 
of C v into the model significantly impacts the recov- 
ery of the linear redshift-space distortion parameter, /3, 
whereas the distance constraints remain unaffected by the 
marginalization over the radiative transfer effects. This 



differs from Wyithe & DijkstraN201l|) where inclusion of 



scale-dependent ionizing background fluctuations lead to 
reduced distance constraints. 

With sufficiently tight priors on C v , the constraints 
on P approach the results given in the previous section, 
as expected. If we have a poor understanding of C v , how- 
ever, the ability to recover j3 drops by an order of mag- 

2 Throughout this work we define the 'perfect knowledge on 
CV as a prior set to o~c = 0.0001, which is to ensure that 
our Fisher matrix elements remain finite, yet still mimic the 
behaviour for perfectly understood parameters. 



nitude. See the left panel of Figure [2] for a graphical rep- 
resentation of the effect of the priors. 

On the other hand, the distance constraints are un- 
affected by the degeneracy between fi and C v . (Compare 
the middle panel of Figure[l]with the left panel.) This is 
because /3 and C v enter into the power spectrum in the 
same way: as far as the distance scales are concerned, it 
makes no difference whether one ma rgin alizes over j3 in 



Equation [2] or /3(1 — C v ) in Equation 
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3.3 LAE galaxy power spectrum 



Wyithe & Dijkstra (20111 show that the magnitude of 



the Lycv radiative transfer parameters varies significantly 
depending on the LAE model considered. In particular 
the magnitude of the effect is significantly larger in the 
absence of a galactic outflow, so that the absorption is 
dominated by infalling IGM. For illustration we consider 
this 'infall' model with an escape fraction of 10 per cent, 
as this is the model with the largest magnitude Lya ra- 
diative transfer effects. Hence, when we include the ra- 
diative transfer effects into our model, we set Cr = 0.05, 
C p = -0.39, and C v = 0.11 in Equation fl3l 
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In Table [3] we provide estimates for the recovery of 
cosmological parameters when we include the radiative 
transfer effects. The linear bias is modified from its fidu- 
cial value of 61 = 2.2 to 61 = 1.9, acting to reduce the 
observed clustering of LAE galaxies. As the level of the 
shot noise is the same, a reduced effective bias implies a 
lower signal-to-noise for measuring the power spectrum 
of LAE galaxies. As a result, the expected constraints 
on the angular diameter distance and the Hubble rate 
are worse than the previous two cases. One can see this 
clearly in the right panel of Figure [T] 

The same is true for /3: due to a smaller effective 
bias, the fractional precision by which we can determine 
f3 is slightly worse than the previous cases. (One can see 
this by comparing the rows of '/3' in Table [2] and [3]) Note 
also that, as the fiducial value of 61 is different due to 
the radiative transfer parameters, the fiducial value of 
/3 = f/bi is also different. 

In Figure [2] we show the one-dimensional likelihood 
distributions for /3 for both the fiducial LAE galaxy model 
(left panel) and the LAE galaxy model (with Lya effects, 
right panel). As already described in the previous section, 
the recovery of /3 from these models is highly sensitive 
to the priors on C v , with marginal improvement on the 
priors breaking the degeneracy between /3 and C v . In the 
right panel, we compare the likelihood distribution for 
the case of perfect knowledge of C v for the fiducial LAE 
galaxy model (the dashed line) to the LAE galaxy model 
(the solid line) , showing the degree of offset that the Lya 
radiative transfer parameters have on the fiducial value 
of /3. 

Thus inclusion of the Lya radiative transfer param- 
eters impacts the recovery of cosmological constraints, 
most notably the growth rate of structure / through the 
recovery of the linear redshift-space distortion parameter 
/3. Unless we have a good prior knowledge on the value of 
C v , it seems hopeless to determine /3 with any precision. 
Fortunately, one can break the degeneracy between j3 and 
C v by including the three-point function (bispectrum), as 
we shall show next. 



4 BISPECTRUM AND NON-LINEAR 
CLUSTERING OF LAE GALAXIES 

If primordial perturbations are Gaussian, linear density 
fields are also Gaussian, in which case the bispectrum of 
linear density fields vanishes. The bispectrum is defined 
as 

(8(k 1 )5{k 2 )8(k s )} = {2Tv) 3 B{k 1 ,k 2 ,k s )S D (k 1 +k 2 +k s ).{19) 

However, non-linear gravitational evolution of den- 
sity fields, and non-linear gravitational and non- 
gravitational evolution of galaxy bias make the observed 
galaxy density fields non-Gaussian. As a result, the ob- 
served bispectrum does not vanish, providing information 
regarding non-linear evolution of density fields. 

In the previous section, we have considered the 
linear-theory power spectrum model for the LAE galaxy 
population. Structure formation is inherently a non-linear 
process, and by considering the true non-linear galaxy 



power spectrum one would expect to increase the con- 
straining power. However, the non-linear power spectrum 
alone will not achieve this, due to the additional parame- 
ters required to fully describe it. Hence the simple linear 
LAE galaxy power spectrum is preferred instead of the 
increased model complexity provided by the non-linear 
LAE galaxy power spectrum. On the other hand addi- 
tional information on large scales may be contained in 
the non-Gaussianity associated with structure formation. 

We use this information to break the degeneracy be- 
tween cosmological parameters and Lya radiative trans- 
fer parameters. There are three effects: (I) gravitational 
evolution of matter density fields, (2) gravitational and 
non-gravitational evolution of galaxy formation (cap- 
tured by galaxy bias), and (3) non-gravitational Lya ra- 
diative transfer effects. 



4.1 Eulerian perturbation theory 

First, we summarize the non-linear gravitational evolu- 
tion of density fields. Specifically, we apply standard Eu- 



lerian perturbation theory (Bernardeau et al. 2002 and 



references within) which, at larger redshifts, has been 
shown to describe the power spectrum measured from 
N-body simulations accurately ( |Jeong fc Komatsu|2006| ) . 

The next-to-leading order corrections to the matter 
density field, S, as well as to the velocity-divergence field, 
7], are generated from the following expressions: 



S(k,z) = X>"( 2 



d 3 9l /'d 3 q n _ 1 



(27T) 



(27T 



/n 
d 3 qJ D (k- 1 £q i )Ft ) (q 1 ,q 2 ,...,q n ) 
i=l 



X5 1 (q 1 )6i{q 2 )...5i{q n ) 



(20) 



' 8 - r d 3 gn _ i 



(2tt) 3 J (2tt) 3 

n 

d 3 q n S D (k - q,)G ( n s) (q 1 ,q 2 , q n 



xS 1 (q 1 )S 1 (q 2 )...8 1 (q n ), 



(21) 



where D(z) is the linear growth factor describing the evo- 
lution of the linear density field, 8i(q i ), which is a Gaus- 
sian random field, and F„ and G„ are symmetrized ker- 



nel expressions generated from recursive relations (Jain 



& Bertschinger 1994 1. We deal only with the next-to- 



leading order expressions, for which the kernels are well 
known: 



F 2 (<h,9 2 )- 7 +7 q2q2 +-3- 



q 2 ) 



3 4 ( gl • g 2 f 
7 7 q\q\ 



1 

¥2 
1 



9 I 9 

qf q 2 



(22) 
.(23) 



4.2 Galaxy bias 

Galaxies are biased tracers of the underlying dark matter 
density field I Kaiser|I 984). Pushing into the weakly non- 
linear regime, we anticipate contributions from both the 
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linear and non-linear mapping of galaxies to the dark 
matter field. 

The bias of galaxies differs from population to pop- 
ulation, and their exact value depends on the underlying 
galaxy formation processes. Typically we expect a scale- 
dependant bias relating the clustering of the galaxies to 
the underlying matter density on small scales, but on 
large scales we expect the bias to be scale-independent. 

To estimate the clustering of galaxies, we Taylor- 
expand the fluctuations in the number density of galaxies, 
8 g (x), in terms of the underlying matter density field 



fluctuations (Fry fc Gaztanaga|1993| |McDonald||2006[ ): 



S g (x) = e(sc) + bi6(x) + -b 2 8(x) 2 + ... 



(24) 



where S(x) is the non-linear matter density field, and 
foi and b 2 are the linear and non-linear bias parame- 
ters, respectively. The e(x) term is a stochasticity pa- 
rameter describing the non-deterministic relationship be- 
tween galaxies and the underlying matter distribution 

JYoshikawa et al. 2001 1. We shall assume that e is a Gaus- 
' _ . . . . ' ' ,0, 



sian field which is not correlated with 8, i.e., (e ) = and 
(e8) = 0. Under this assumption, e does not contribute 
to the bispectrum, and thus we shall ignore stochasticity 
throughout this paper. 



4.3 LAE kernel expressions 



In Section [2T2| we outline the derivation for the first-order 
Lya radiative transfer effects. To derive higher-order ex- 
pressions for the Lya radiative transfer effects, we Taylor- 
expand nL yQ1 (> Fq) about the three non-gravitational 
Lya radiative transfer effects, in analogy to the galaxy 



bias derivation (Fry & Gaztanaga 1993). We obtain 



n Lya {> L , po, To, 5{x)) = <n$ a (l + n$ a + n {2) ya 



(25) 



where n 



(i) 



and n 



(2) 



are the first- and second- 



'Lya "Lya 

order Taylor-expanded expressions, respectively, evalu- 
ated around the mean quantity n^ a . The term n L y a is 
given in Equation [fi] and n^ya is given by 



-(2) 



1 1 



2n (0) 



(r-r ) 



or 



+ (P~ Po) 



+ 



d(ar c 



H 



F ,T 

d 



dp 



Fq,P0 



a 



d(ar com ) 



(26) 



The number density of LAE galaxies above a flux limit 
can then be written as 



n Lya (> F )^n^ a [l + 8 g (x)} [ 



where we have substituted Equation [25] into Equation [4] 
To proceed further, we firstly expand the expressions 
in Equations [6] and |26[ and then recast the above deriva- 
tives as explicit constants with respect to their radia- 
tive transfer effect. Once this has been performed, we 
can rewrite Equations [fj] and [26] as 



n (1) 



= 5 T (x)Cr + 8 p (x)C p + 6 v (x)C v 



(28) 



and, 



"Lya 



- [CrrSl(x) + Cp P S 2 p (x) + C vv 8*(x) 

+Cr P Sr(x)5p(x) + Crv8r(x)S v (x) 
+C pv 8 p (x)S v (x). 



(29) 



Equations [28] and [29] contain the first- and second- 
order Lya radiative transfer coefficients. In Appendix [A] 
we derive the explicit expressions for the first-order Lya 
radiative transfer coefficients, and use the same basic 
ideas to also derive the second-order coefficients. 

To the second order, the fluctuations in the number 
density of LAE galaxies relative to the mean are given by 



8 Lya (x) = [l + 8 g (x)] [ 



1 + n (1) +n (2) 1 

1 ~ "''Lya T "Lya 



(30) 



To generate the Lya radiative transfer kernels that de- 
scribe the modification to the galaxy power spectrum and 
the higher-order corrections, we expand Equation |30| up 
to the second order in fluctuations, and take the Fourier 
transform of the corresponding expression ( Appendix |B|. 
In the resulting expression, the fluctuations of the local 



density, 8 p (x), are given by Equation 20 The fluctuations 
in the line-of-sight peculiar velocity field in Fourier space 
can be expressed as 8 v (k) = — ffi 2 r/(k), where r](k) is the 
fluctuation in the velocity field given by Equation |21| The 
expression due to fluctuations in the UV background can 
be somewhat more complicated and is outlined in Ap- 
pendix [B] 

In Appendix |B] we also derive the redshift-space ex- 
pressions for the Lya radiative transfer effects. Finally, 
in analogy with Equations |20| and |21| we write the fluc- 
tuations in the number density of LAE galaxies in real 
space as 



.CM) = trwffef*^ 

n 

d a q n 8 D {k-Y J <lr)Z { n S) {q 1 ,q 2 



3Ly 



x8i{q 1 )8 1 {q 2 )...8i(q n 
and those in redshift-space as 

-13 



(31) 



ra=l J V 



d 3 9l f<Tq 



(2tt) 3 

X f d 3 q n 8 D (k-^ qi )K^( qi , q2 , 

J i=l 

x5i(q 1 )6i(q 2 )...5i(q n ), 



(32) 



where the Lya radiative transfer kernels, and K„\ 
are given in Appendix [B] 



4.4 LAE bispectrum and reduced bispectrum 

The second-order terms in Equation [32] yield a non- 
vanishing bispectrum of the fluctuations in the number 
density of LAE galaxies in redshift space: 

B hya , s (ki,k 2 ,k 3 ) = 2\K[ s) (k 1 )K[ s \k 2 )K^\k l ,k 2 ) 
xP L (k 1 )P L (k 2 ) + (2 eye.)], (33) 
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where K[ s) (k) and i^ s; (fci, &2) are found in Appendix 

m 

The bispectrum in redshift space depends on six vari- 
ables: three wavenumbers, ki, &2, and ks, giving the sides 
of a triangle; and the cosines of the angles that these three 
vectors make with the line-of-sight direction, (ii, (12, and 
(j,3. However, due to the triangular condition, these six 
variables are not all independent. Instead, the bispectrum 
can be written as a function of five independent variables 
(Scoccimarro et al. 1999 Smith et al. 20081: three pa- 



rameters (fei, k2, and the angle between them, cos(#i2)) 
define the shape of the triangle; and the remaining two 
parameters (/ii and (j>) define the orientation of the tri- 
angles with respect to the line-of-sight. 

By convention, we align the first wavevector, ki, to 
the line-of-sight direction, z, about which the triangle 
can be rotated through the azimuthal direction ((f)). The 
cosines of the angles that k^ and £3 make with the line- 
of-sight direction are given by 



fi2 



M3 = 



^iicos( 

fcl 



/i^sin(6'i2)cos(0) 



k. 



(34) 
(35) 



Here the last equality comes from the triangular condi- 
tion, fcl + k2 + &3 = 0. 

The 'reduced' bispectrum is given by the ratio of the 
bispectrum to the products of the power spectra: 



<2Lya,s(fcl, fe2, fe3 



i?Ly a ,s(fcl, &2, fc3 



PLyc«,s(fcl)-PLya, S (fc2) + 2 CyC. 



(36) 



This quantity is insensitive to the overall amplitude of 
the power spectrum, as the second-order expression for 
the bispectrum given in Equation [33] is proportional to 
the products of the power spectra. This properly removes 
the degeneracy between the galaxy bias parameters and 
the amplitude of the matter power spectrum. 

4.5 Fisher matrix 

Before generating the expected cosmological constraints 
using the Fisher matrix, let us first summarize the model 
parameters characterizing the higher-order (non-linear) 
terms. 

It is important to note that, unlike the previous work 
which simply multiplies the real-space bispectrum by the 
linear redshift distortion factors ( Scoccimarro et al.|1999 



Sefusatti et al.||2006| [Sefusatti fc Komatsu||2007| ), we in- 
clude the full wavenumber dependence of the redshift- 
space distortion up to the second order. By including the 
full second-order redshift-space distortion, we gain addi- 
tional information which helps to further break the de- 
generacies between the cosmological information and the 
radiative transfer effects, especially those associated with 
the velocity gradient. 

The second-order redshift-space Lya kernel, after re- 
moving the scale dependence of the ionizing background 
effect, is given by (see Appendix IB} 



^ s) (fci,fc 2 



1- 



~&2- ^/( M i+M2)C + 6 1 F 2 (s) (fc 1 ,fc 2 ) 



k2z 
fe| 



+ ^(fcl2Ml2./) 2 



r [h - ft-£c v ] + 

kizk 2z 



£-2 £.2 
^1^2 



- a)G 2 s) (fc!, fca) + \fnl£c vv 

+ ^6l(fcl2/W) 



ki z 

fc? 



+ 5 (fcia/iia/) a (l-a) 



k\ 



kizk 2 z 



2i.2 
2 



(37) 



where we define 



C 

kij 
ki z 



bi (Crr + 2C r ) + C Pi 
biC v + C pv + biCrv, 

| ki ~\~ kj | j 

( k i -J - kj ) • Z 

| ki -\- kj | 
ki • z. 



&2(l + C r ) + 26 1 (Cr P + C p ), 
(38) 
(39) 
(40) 

(41) 

(42) 



Here, 62 is the effective non-linear galaxy bias modified by 
various first- and second-order radiative transfer effects, 
and fei is defined in Equation 
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In our model, we choose to keep the second-order 
effect due to the peculiar velocity gradient (C vv ) sepa- 
rate, as this has the potential to be degenerate with the 
growth rate of structure, /. Additionally, C, which con- 
tains the linear-order effects with respect to the peculiar 
velocity gradient, can also become degenerate with /. We 
note that setting the Lya radiative transfer effects to zero 
reduces Equation [37] to the typical second-order redshift- 
space galaxy kernel, as required. 

To generate the expected constraints on the cosmo- 
logical parameters, we calculate the Fisher matrix for 
both the bispectrum and the reduced bispectrum. The 
Fisher matrix for the bispectrum is 



1 dB g (k,n)dB g (k,iJ,) 



80, 



f « - E 

ki,k 2 ,k 3 <k r) 

and for the reduced bispectrum is 

1 dQ g (k,fj.) dQ g (k,n) 

c 

ki,k2,k 3 <k ma:r 



Fii = E 



(43) 



(44) 



where the Fisher matrices are summed over all possible 
triangular configurations. 

The variances for the bispectrum and the reduced 
bispectrum are given, respectively, by 



2 SsVsurvey 

ob = 



and 



2 SsVsurvey 

a Q = 



Ptat(kl)P t „t(k2)Pt t(k 3 ) 

Ptot(kl)Ptot(k 2 )Ptot(k3) 



(45) 



(46) 



N t [PLy Q , S (fcl)P L y Q , S (fc2) + 2cyC.] 2 ' 

Here Sb is the symmetric factor describing symmetry of 



+// i i2(l C V )G 2 (fci,fc2)+ / /I1P2C1™ the side lengths of a given bispectrum triangle (sb = 
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6, 2, 1 for equilateral, isosceles, and general triangles, re- 
spectively) and P to t (k) is the sum of the power spectrum 
component and the Poisson shot noise: 



Ptot{k) = P Ly «, s (fc) + — , 



(47) 



where n g is the number density of LAE galaxies. The 
quantity N t is the total number of available triangles: 

V B 



N t = 



(48) 



where &f is the fundamental frequency and 

if hi / kj + kk , 



V B = 2nd^d(j>k 1 k2k3(Aky i x 



if ki — kj -\- kk ■ 



(49) 



This states that, for 'collapsed' (or 'co-linear') triangles 
defined by ki = kj + kk, the bispectrum volume is re- 
duced by a factor of two. (For the derivation of Vb, see 
Appendix [CJ. In the simplest spherically averaged sce- 
nario this reduces to 



V B = 32ir 2 k 1 k 2 k 3 (Ak) 3 x 



if ki ^ kj -\- kk i 
if ki — kj -\- kk ■ 



(50) 



5 CONSTRAINTS FROM THE REDUCED 
BISPECTRUM ALONE 

In Section [3] we show that the growth rate of structure, 
/ (or /3), is completely degenerate with the Lya radiative 
transfer effect due to the velocity gradient, C v , as long as 
we rely only on the power spectrum. 

However, the bispectrum provides additional con- 
straining power that can be used to break the degen- 
eracy between the growth rate of structure, /, the lin- 
ear galaxy bias, b\, and C v . This is because, unlike the 
power spectrum which just tells us the amplitude of the 
fluctuations at a given scale, the bispectrum tells us also 
how the structure forms. For example, one needs infor- 
mation in the bispectrum in order to reproduce the 'cos- 
mic web,' the filamentary structures in the universe. The 
power spectrum cannot distinguish between the distribu- 
tion with random phases and that with the filamentary 
structures, as it is sensitive only to the amplitude of the 
fluctuations. As a result, the bispectrum can distinguish 
between the structures caused by gravitational and non- 
gravitational effects. 

In this section we firstly generate the expected cos- 
mological constraints from the reduced bispectrum. As 
mentioned previously, the reduced bispectrum is insen- 
sitive to the amplitude of the matter power spectrum. 
We again consider the same two models; a fiducial model 
where we set the radiative transfer coefficients to be zero 
but marginalize over them, and a model where we use ex- 
plicit values for the Lya radiative transfer effects in our 
redshift-space expressions. 

Since the recovery of the growth rate of structure / 
is most affected by the radiative transfer effects, we inves- 
tigate the two-dimensional joint likelihood distributions 
for / with each of the other model parameters (marginal- 
ized over all the remaining model parameters). For the 
remainder of this work, we set the non-linear galaxy bias 
to be b 2 = 1.5. 



5.1 Fiducial LAE reduced bispectrum 

We first consider our fiducial model where we set all Lya 
radiative transfer coefficients to zero, but marginalize 
over the Lya effects. With the addition of the bispectrum, 
the number of parameters in our model has increased to 
eight. These include three cosmological parameters: /, 
Iii(Da), and \n(H); three radiative transfer parameters: 
C v , C vv , and C; and the linear and non-linear galaxy 
biases: b\ and b 2 . 

We find that the constraints generated from the 
reduced bispectrum contain no strong degeneracies be- 
tween / and the radiative transfer parameters (see Fig- 
ure [3]). The reduced bispectrum does however exhibit 
some degeneracies between /, and the galaxy bias pa- 
rameters, 61 and b 2 . 

To understand this result, let us write the reduced 
bispectrum given in Equation |36| as 

QLya,s(fel, ^2, k^) 

2K[ S) (fci) K[ s) (k 2 )Ki s) (fci ,k 2 )P L (fci)Pi (fe) + 2 eye. 



[k[ 3) {k 1 )Y[K\ s, (k 2 )YP L {k 1 )P L {k 2 ) + 2cyc. 



where 



^ S) (fcl,fe2 



^-K[ s) (k) = l + p fJ 2 (l-C v 
bi 



(51) 
(52) 



^ s) (fci,fc 2 



1^ - , l M + A)c + ^ (s) (fc 1 ,fe2) 
+/3 M ?2(i-a,)G< s) (fc 1 ,fc 2 )] 

CO 4 )- (53) 



+ ^/3(fcl2A*12) 



k\z k 2z 

fv 1 fi, 2 



Here, C(/i 4 ) are the term s that contain four powers of 
cosines (see Equation 37 for the full expression of K^). 
These terms contribute less, as their contributions are 
important only near the line-of-sight direction, for which 
the number of available modes is limited. 

Equations |51[ |52| and [53] show that the reduced bis- 
pectrum determines the following parameter combina- 
tions: 

• 61 from the overall amplitude of the first four terms 
in K. 



b 2 /b\ from a constant, fc-independent term in K. 



/3(1 — Cv) from and the term proportional to 



• PC from the second term in and 

• /3 from the last term before 0(fj, 4 ) in K% . 

Recalling j3 = f/bi, there are five unknown variables (61, 
b 2 , f, C v , and C), and the reduced bispectrum yields five 
combinations of these variables. 

From the Fisher matrix calculations, we find that 
the reduced bispectrum primarily yields b 2 /bi and /3. The 
information on 61 coming from the first four terms in A" 2 S ' 
breaks a complete degeneracy between b 2 and &i and /, 
but correlations between these parameters still remain. 
One can see this in the first two panels in Figure [3] On 
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/ / / 

Figure 3. Two-dimensional joint marginalized likelihood distributions computed from the fiducial LAE galaxy reduced bispectrum 
alone (no Lya radiative transfer effects, but including marginalization over C v , C vv , and C). We show the correlations between 
the growth rate of structure, /, and various parameters including (clockwise from top left): the linear bias, hi, non-linear bias, 02, 
linear peculiar velocity Lya effect, C v , non-linear peculiar velocity Lya effect, C vv , the non-linear combination of other radiative 
transfer effect, C, angular diameter distance, ln(D^), and the Hubble rate ln(H). The solid and dashed curves show the 1- and 
2-cr joint marginalized constraints, respectively. 




0.8 1 1.2 0.8 1 1.2 0.8 1 1.2 

/ / / 

Figure 4. One-dimensional marginalized likelihood distributions for the growth rate of structure, /, for the fiducial case (no 
Lya radiative effects added, but including marginalization over C v , C' vv and C) generated from; {left) the LAE galaxy reduced 
bispectrum only, {centre) the LAE galaxy power spectrum combined with the LAE galaxy bispectrum, and {right) the LAE galaxy 
power spectrum combined with the LAE galaxy reduced bispectrum. The various curves denote different priors added to C v \ black 
solid: Perfect knowledge of C' v , grey dotted: o~q = 0.01, grey dot-dashed: cfq^ = 0.1, grey dashed: gq^ = 0.5, and grey solid: no 
priors added. 



the other hand, we do not find much correlation between 
/ and the radiative transfer parameters, C, C v , and C vv 
(see the third to fifth panels of Figure [3|. 

While the reduced bispectrum does break the degen- 
eracy between / and C v seen in our power spectrum anal- 
ysis, it cannot provide a strong constraint on /. In the 
fourth column of Table [4] we provide the 1-a constraints 
generated from the one-dimensional likelihood distribu- 
tion for /. In the fiducial case with no additional priors, 
we find the 1-a constraint on the growth rate of structure, 
/, to be 0.16 (17 per cent). It is important to note that, 
while the constraints are relatively weak, they are on / 
as opposed to /5. Also, /3 and C v are totally degenerate 
in the LAE power spectrum, and thus the error bar on 
/3 is infinite unless we put a prior on C v . Therefore, the 
reduced bispectrum provides a massive improvement on 



the constraint on /: the error bar shrinks from infinity to 
17 per cent. 

In the left panel of Figure |4j we show the one- 
dimensional likelihood distributions for the growth rate 
of structure, /, for various priors on C v . The addition of 
priors to C v does not improve the constraints on / from 
the reduced bispectrum alone, as the reduced bispectrum 
contains no degeneracy between / and C v . 

In the fourth column of Table [4] we also provide 
the 1-a constraints from the one dimensional likelihoods 
for Iii(Da) and ln{H) given various priors on C v . We 
find that (independent of priors on C v ) the fiducial LAE 
galaxy reduced bispectrum can recover the angular diam- 
eter distance scale at 3 per cent and the Hubble rate at 2.5 
per cent. This should be contrasted with the 1.1 per cent 
and 1.3 per cent errors on Da and H expected from the 
fiducial LAE galaxy power spectrum. Clearly the reduced 
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Priors on C v 


Parameter 


Model 


R BS 


PS + BS 


PS + R BS 








1-(t (per cent) 


Tcr (per cent) 


1-<t (per cent) 


1\U 1 1 1U1 


J 


T-Ti nn^i £i 1 
_L 1U LlLlctl 


1 fid5 ('\ (\ 
u.ru'itj ^iu.y^^ 






Perfect knowledge 


I 


Fiducial 


0.1635 (16.82) 


0.0507 (5.21) 


0.0565 (5.81) 


0.01 


f 


Fiducial 


0.1635 (16.82) 


0.0520 (5.35) 


0.0576 (5.93) 


0.1 


f 


Fiducial 


0.1636 (16.83) 


0.1055 (10.85) 


0.1063 (10.93) 


0.5 


f 


Fiducial 


0.1642 (16.89) 


0.1555 (16.00) 


0.1535 (15.79) 


No Priors 


ln(D A ) 


Fiducial 


0.0303 (3.03) 


0.0076 (0.76) 


0.0103 (1.03) 


0.01 


ln(D A ) 


Fiducial 


0.0301 (3.01) 


0.0075 (0.75) 


0.0101 (1.01) 


0.1 


hi(D A ) 


Fiducial 


0.0302 (3.02) 


0.0075 (0.75) 


0.0102 (1.02) 


No Priors 


ln(H) 


Fiducial 


0.0251 (2.51) 


0.0084 (0.84) 


0.0115 (1.15) 


0.01 


ln(H) 


Fiducial 


0.0249 (2.49) 


0.0083 (0.83) 


0.0112 (1.12) 


0.1 


HH) 


Fiducial 


0.0249 (2.49) 


0.0083 (0.83) 


0.0113 (1.13) 



Table 4. We show the 1-<t constraints expected from the reduced bispectrum (R BS), the power spectrum combined with the 
bispectrum (PS + BS), and the power spectrum combined with the reduced bispectrum (PS + R BS). No Lya radiative transfer 
effects are included, but the likelihood is marginalized over C v , C vv , and C. The first five rows show the l-cr constraints on / for 
various priors on C v , after marginalizing over 61, 62, C, C' v , C vv , ln(D^), ln(H), and the amplitude [ln(j4)]. The last six rows 
show the l-cr constraints on the distance parameters, ln(D A ) and ln(H), marginalized over the remaining model parameters. 



Priors on C„ 


Parameter 


Model 


R BS 
l-<7 (per cent) 


PS + BS 
l-cr (per cent) 


PS + R BS 
l-cr (per cent) 


No Priors 

Perfect knowledge 

0.01 

0.1 

0.5 


/ 
/ 
/ 
/ 
/ 


LAE effects included 
LAE effects included 
LAE effects included 
LAE effects included 
LAE effects included 


0.2104 (21.64) 
0.2089 (21.49) 
0.2089 (21.49) 
0.2090 (21.50) 
0.2099 (21.59) 


0.2039 (20.97) 
0.0632 (6.50) 
0.0645 (6.64) 
0.1262 (12.98) 
0.1965 (20.21) 


0.2014 (20.72) 
0.0672 (6.91) 
0.0685 (7.05) 
0.1270 (13.06) 
0.1943 (19.99) 


No Priors 

0.01 

0.1 


ln(D A ) 
ln(D A ) 
ln(D A ) 


LAE effects included 
LAE effects included 
LAE effects included 


0.0378 (3.78) 
0.0376 (3.76) 
0.0376 (3.76) 


0.0099 (0.99) 
0.0098 (0.98) 
0.0099 (0.99) 


0.0120 (1.20) 
0.0119 (1.19) 
0.0119 (1.19) 


No Priors 

0.01 

0.1 


ln(H) 
ln(H) 
ln(H) 


LAE effects included 
LAE effects included 
LAE effects included 


0.0305 (3.05) 
0.0302 (3.02) 
0.0303 (3.03) 


0.0107 (1.07) 
0.0106 (1.06) 
0.0107 (1.07) 


0.0133 (1.33) 
0.0131 (1.31) 
0.0132 (1.32) 



Table 5. Same as Table[4] but for the first-order Lya radiative transfer effects given by Cr = 0.05, C p = —0.39, and C v = 0.11. 



bispectrum alone provides weaker distance constraints. 
This is not surprising, as the distance information is con- 
tained in the shape of the power spectrum (e.g., baryon 
acoustic oscillation (BAO) and Alcock-Paczynski (AP) 
test), which is largely divided out in the reduced bispec- 
trum. 



the precision with which we can constrain the growth rate 
of structure / has been reduced to an error of 0.21 (22 
per cent) compared to 0.16 (17 per cent) for the fiducial 
model. Once again, this is due to the reduced effective lin- 
ear galaxy bias, which reduces the signal-to-noise ratio of 
the LAE power spectrum relative to the shot noise. 



5.2 LAE reduced bispectrum 

We now consider the inclusion of Lya radiative transfer 
effects by adding the linear Lya radiative transfer model 
parameters, Cr = 0.05, C p 



-0.39, and C v 



0.11 



from Wyithe & Dijkstra (20111. The inclusion of these 
parameters modifies the effective bias parameters, 61 and 
&2, and the Lya radiative transfer effects associated with 
C. We still set the fiducial values of the second-order Lya 
radiative transfer coefficients to vanish. Although we set 
Cvv = 0, we still marginalize ov6r Cvv 

in our models. 
In the fourth column of Table [5j the l-cr constraints 
on /, \n(D A ), and ln(//) are generated from the likeli- 
hood distributions for various priors on C v as per the 
previous section. With the inclusion of the Lya effects, 



6 COSMOLOGICAL CONSTRAINTS FROM 
COMBINING THE POWER SPECTRUM 
AND BISPECTRUM 

We next discuss the improvements on the cosmological 
constraints available when we combine the LAE power 
spectrum with either the LAE reduced bispectrum or 
the bispectrum. When combining the reduced bispectrum 
(and the bispectrum) to the information from the power 
spectrum, we assume that there is no covariance between 
the power spectrum and the reduced bispectrum (or the 
bispectrum), which is incorrect. Therefore, the numeri- 
cal values of the l-cr constraints on various parameters 
reported here should be considered as lower bounds. 
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Figure 5. Two-dimensional joint marginalized likelihood distributions computed from the fiducial LAE galaxy power spectrum 
combined with the fiducial LAE galaxy reduced bispectrum (no Lya radiative transfer effects added, but including marginalization 
over C v , C vv , and C). We show the correlations between the growth rate of structure, /, and various parameters including 
(clockwise from top left): the amplitude, ln(A), linear bias, bi, non-linear bias, 621 linear peculiar velocity Lya effect, C v , non- 
linear peculiar velocity Lya effect, C vv , the non-linear combination of other radiative transfer effect, C, angular diameter distance, 
ln(D^), and the Hubble rate, ln(if). The solid and dashed curves show the 1- and 2-a joint marginalized constraints, respectively. 



6.1 Fiducial LAE model for power spectrum 
and bispectrum 

We first consider our fiducial model in which all Lya ra- 
diative transfer coefficients are set to zero. The number of 
parameters in this model is nine. While the LAE galaxy 
reduced bispectrum is insensitive to the amplitude of the 
matter power spectrum, we must marginalize over the 
amplitude information in the LAE galaxy power spec- 
trum. The parameters include four cosmological param- 
eters: the amplitude [1h(j4)], /, \n(T>A), and ln(H); three 
radiative transfer parameters: C v , C vv , and C; and the 
linear and non-linear galaxy biases: b\ and 62. 

6.1.1 Combined power spectrum and reduced 
bispectrum 

Figure[5]shows the expected constraints from a joint anal- 
ysis of the reduced bispectrum and the power spectrum 
on various pairs of parameters involving /. Comparing 
this figure with Figure [3] we find that adding the power 
spectrum does not improve the constraints on / and the 
bias parameters very much, but improves the constraints 
on all the other parameters. Figure [6] shows this more 
clearly: adding the power spectrum information does not 
improve the constraint on /, but it substantially improves 
the constraint on C v . 

What does this imply? This implies that the uncer- 
tainty in / is now dominated by the correlation between 
/ and the bias parameters - the correlation that we have 
discussed in Section [5.1| Comparing the fourth and sixth 
columns of Table [4] shows this quantitatively. 

Comparing the fourth and sixth columns of Table [4] 
also shows that adding the power spectrum does improve 
the constraints on Da and H substantially, as the power 
spectrum contains features such as BAO and AP test, 
whereas such information is largely cancelled out in the 
reduced bispectrum. 




Figure 6. Comparison of the joint two-dimensional con- 
straints on / and C v . Outer two ellipses correspond to the 1- 
and 2-(T constraints generated from the fiducial LAE galaxy 
reduced bispectrum only. Two narrower ellipses correspond 
to the 1- and 2-cr constraints generated from the fiducial LAE 
galaxy power spectrum combined with the fiducial LAE galaxy 
reduced bispectrum. 



Nevertheless, as the reduced bispectrum still has 
some sensitivity to these features (i.e., cancellation is not 
exact), the constraints on Da and H from the power 
spectrum and the reduced bispectrum are slightly better 
than those from the power spectrum alone. Comparing 
the third column of Tableland the sixth columns of Ta- 
ble[4] we find that the expected constraints improve from 
1.1 to 1.0 per cent for Da and 1.3 to 1.2 per cent for H. 



6.1.2 Combined power spectrum and bispectrum 

Next, we combine the power spectrum with the bispec- 
trum (rather than the reduced bispectrum). As far as / 
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Figure 7. Two-dimensional joint marginalized likelihood distributions for the angular diameter distance, Da, and the Hubble 
rate, H. Shown also are the 1- (solid) and 2-a (dashed) joint likelihood contours. From left to right show the constraints generated 
from the LAE galaxy reduced bispectrum only, the LAE galaxy power spectrum and bispectrum combined, and the LAE galaxy 
power spectrum and reduced bispectrum combined. Top panels: Fiducial case, with no Lye? radiative transfer effects added to the 
fiducial parameters, but marginalizing over C v , C vv , and C. Bottom panels: The inclusion of the first-order hya radiative transfer 
effects, C v = 0.11, Cr = 0.05, and C p = —0.39, and marginalizing over Cv j Cvv ■) <uid C . 



is concerned, we have the same story: adding the power 
spectrum does not improve the expected error bar on / 
(see the fifth column of Table HJ. 

On the other hand, a joint analysis of the power spec- 
trum and the bispectrum yields a significant improve- 
ment on the angular diameter distance and the Hubble 
rate. This is because the bispectrum also contains the 
BAO features and the AP test in its wavenumber de- 
pendence. However, this could be due to our ignoring a 
covariance between the power spectrum and the bispec- 
trum: a correlation between them would degrade the con- 
straints in a joint analysis. This point requires a further 
investigation. 



6.2 LAE model for power spectrum combined 
with bispectrum 

Finally, we consider the inclusion of Lya radiative trans- 
fer effects on the recovery on /, ln(D^), and \xi(H), 
by adding the linear Lya radiative transfer parameters, 
Cr = 0.05, C p = -0.39 and C v = 0.11 from |Wyithe~fc| 
Dijkstra (2011). Table [5] shows the results: the expected 
constraints are slightly weaker than those from the fidu- 
cial case, which is again due to a smaller effective bias, 
6i, reducing the amplitude of the signal relative to the 
shot noise. 



7 SUMMARY AND CONCLUSION 

In this paper we have studied how the radiative transfer 
effects alter the power spectrum and bispectrum of LAE 
galaxies, and how we can use these properties to separate 
the radiative transfer effects and the cosmological effects, 
so that we can improve the cosmological constraints de- 
rived from them. 



First, as a follow up to Wyithe & Dijkstra ( 2011 1, we 



show that the growth rate of structure (/) and the param- 
eter C„ describing the radiative transfer effects of velocity 
gradients are completely degenerate in the linear power 
spectrum. Next, by performing a perturbation theory ex- 
pansion of the Lya radiative transfer effects, we derive 
the next-to-leading order corrections to the density fields 
of LAE galaxies. This allows us to derive the leading- 
order expression for the bispectrum of LAE galaxies. We 
then show that the reduced bispectrum alone can deter- 
mine / and C v separately, leaving no degeneracy between 
them. Adding the power spectrum information to the re- 
duced bispectrum does not improve the precision of / 
further, as the precision of / is now limited by remaining 
correlations between / and the galaxy bias parameters, 
bi and 62- 

We find that HETDEX-like surveys of LAE galaxies 
can determine / to about 20 per cent accuracy, if we do 
not assume any prior information on C v . Including the 
prior on C„, the uncertainty on / can be reduced down to 
7 per cent. Note that this is the uncertainty on /, rather 
than on /3 = f/bi. 

We find that the constraints on the angular diameter 
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distance and the Hubble expansion rate are not directly 
affected by the radiative transfer parameters (with the 
caveat that we have assumed that the effect of the UV 
ionizing background fluctuation is not scale dependant). 
The only indirect effect is a slight reduction of the effec- 
tive linear galaxy bias, which reduces the amplitude of 
the LAE power spectrum with respect to the shot noise, 
thus slightly increasing the uncertainties in the angular 
diameter distance and the Hubble rate. Comparison be- 
tween the top and bottom panels of Figure [7] shows this 
graphically. 

Finally, to summarize the results of this work, we 
provide Table [6] detailing the constraints on /3, /, Da, 
and H expected from HETDEX-like surveys. This table 
shows how powerful such surveys are in terms of mea- 
suring the distance, the expansion rate, as well as the 
growth rate of the structure in a high-redshift universe, 
and the determination of these quantities are not signifi- 
cantly compromised by the Lya radiative transfer effect. 
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Parameter 


Model 
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PS + BS 
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0.0520 
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LAE effects 
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0.2104 
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LAE effects 
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0.2089 


21.5 


0.0645 
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0.1 


f 


LAE effects 


0.0633 0) 
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0.2090 


21.5 


0.1262 


13.0 


0.1270 


13.1 


No Priors 


\n(D A ) 


Fiducial 


0.0110 


1.10 


0.0303 


3.03 


0.0076 


0.76 


0.0103 


1.03 


No Priors 


\n(D A ) 


LAE effects 


0.0128 


1.28 


0.0378 


3.78 


0.0099 


0.99 


0.0120 


1.20 


No Priors 


ln(ff) 


Fiducial 


0.0132 


1.32 


0.0251 


2.51 


0.0084 


0.84 


0.0115 


1.15 


No Priors 


\n(H) 


LAE effects 


0.0151 


1.51 


0.0305 


3.05 


0.0107 


1.07 


0.0133 


1.33 



Table 6. Summary table of the relevant l-cr constraints (and fractional errors) generated from the one-dimensional likelihood 
distributions for /3, /, ln(D A ), and ln(A). Models considered are 'Galaxy' (no Lye* radiative transfer effects, and no marginalization 
over the radiative transfer parameters), 'Fiducial' (no Lya radiative transfer effects, but marginalized over the radiative transfer 
parameters), and 'LAE effects' (non-zero Lya radiative transfer parameters are included and marginalize over). We consider a 
power spectrum only model (PS), reduced bispectrum only (R BS), power spectrum combined with the bispectrum (PS + BS), 
and the power spectrum combined with the reduced bispectrum (PS + R BS). We explore three priors on the Lya radiative 
parameter associated with the peculiar velocity, C v : no prior at all, prior of 0.01, and prior of 0.1. 
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APPENDIX A: Lya RADIATIVE TRANSFER 
COEFFICIENTS 

Throughout this work, we denote the Lya radiative trans- 
fer effects as constants, which encompass the derivatives 
of the transmission function with respect to the Lya ra- 
diative transfer effect. Here we outline the derivations for 
the first-order Lya radiative transfer coefficients, from 
which the second-order constants can be easily calcu- 
lated. In the following derivations, each of the Lya ra- 
diative transfer effects is expressed as a function of an 



arbitrary transmission function, T (Wyithe & Dijkstra 
20iT) >. 

First, let us consider the fluctuations in the ioniz- 
ing background, which are taken from the expression in 
Equation [6] 



(r — To) (9n Ly c 



r-r 



or 

FoXo 

£>log(T) 



aiog(r) 



dn 



Lya 



rr ,r 



91og(T) 



F ,To 



SrCr, 



(Al) 



where we have defined 8r 



and 



Cv 



i aiog(r) 



(o) aiog(r) 



'Lya 



dn^yc 



rr ,r, 



dlog(T) 



(A2) 



Fo,T 



Following the same idea, we redefine the other two 
terms in Equation [6] as fluctuations in the density field, 
and in the line-of-sight velocity field. For the density field, 
we find 



n (0) 



(p~Po) 



dnLy 



SpCp, 



(A3) 



Fa, Pa 
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where 



Bl Real-space Lya kernels 



C p = 



To,po 



dlog(T) 



(A4) 



Fo,T 



However for the velocity gradient, we rewrite the deriva- 
tive as 



dv z 



n (0) V 


d(ar com ) 


1 1 


/ dv z 


n (0) a 


Vd(r com ) 


dv pcc 


dn Lya 


d(r CO m) 





H 



driLya 



Ha 



Fo,PO 



3 ( ^ 

U 1 d(ar com ) 

dn^ya 



-Fo.PO 



Fo,PO 



(A5) 



where we have expressed the total velocity as v — 
H(ar com ) + v pcc , yielding the velocity gradient 



dv dVp ec 

— Ha + 



d(r CO m) 



d(r c 



(A6) 



Following the case for the ionizing background, we rewrite 
the partial derivative as 



dUhya 



Qdv 

dr C{ 



Fq,pq 

aiogcn 



aiog(^) 



91og( 



dv z 



Using 

dv z 
d?"coii 

we find 



Tq,po 



Ha, 

1 

Wa 



aiog(T) 



Fo,T 



(A7) 

(A8) 
(A9) 



dv z 
\d(ar coln ) 

Lya 

1 1 d?w 



d dv > 



d(ar com ) 

Slog(T) 



n (o) a Ha d(r com ) dlog(^) 



*Fo,PO 

<9n Lyc 



To,po 



aiog(r) 



where (L = ^ Vpcc and 

u iia a(r com J ' 



Cv 



d\o g (D 



dn^ya 



4t^(fc) To , po ^g(r) Foir ; 



F ,T 

(A10) 



(All) 



APPENDIX B: Lya RADIATIVE TRANSFER 
KERNELS 

Here we outline the derivation of the higher-order real- 
and redshift-space kernels for the fluctuations in the num- 
ber density of LAEs. 



In this subsection, we derive the symmetrized kernels, 
Zn^ , which give the LAE galaxies in real space as 

d 3 9l f d 3 g n _! 
WJ (2tt)3 



5h ya {k, z) 



d S q n S D (k -~^2q i )Zi s) (q 1 ,q 2 , ...,q n ) 



xS 1 (q 1 )S 1 (q 2 )...5 1 (Qn)- 



(Bl) 



Expanding out Equation [30] and keeping terms only 
to second order in fluctuations, 



SLyctix) = S g (x) + n^ a + ^Lya + S 9 ( X ) n hy 



S (S) 



(B2) 

Now, substituting in Equations |24| [28] and |29| into Equa- 
tion |B2| and expanding, again only keeping terms up to 
the second order in fluctuations, 

5hya{x) 

1 

= bi$(x) + ^b 2 S(x) z + S r {x)C r + S p (x)C p + S v (x)C v 

+ \ [CrrSUx) + C pp 5 2 p (x) + C vv S 2 v {x)] 

+Ct p St(x)S p (x) + CrvSv(x)S v (x) + C pv S p (x)5 v (x) 

+bi6(x) [S r (x)C T + 8„{x)G p + S v (x)C v ] . (B3) 

To model the effect of the fluctuating ionizing back- 
ground, we convolve the overdensity of sources (in this 
case the galaxies), 5 g (x), with a f unction of the form 
oc exp[-(cc - x )/\]/[(x - x ) 2 ] (|Morales fc Wyithe 



2010 1, where A is the mean free path of the ionizing pho- 
c find the Foi 

arctan(|fc| A) 



tons. Hence we find the Fourier transform of the ionizing 
fluctuations as 



5 r (k) = S g (k)- 



felA 



which, to the second order, yields 
Sr(k) = 



M(fc) + ;j,M(fc) 2 



arctan(|fc| A) 
IfclA ' 



(B4) 



(B5) 



Taking the Fourier transform of Equation |B3| and 
using Equation |20| |21| and |B5| for the density field, 
peculiar-velocity field, and the ionizing background, re- 
spectively, and again keeping terms up to the second or- 
der only, one finally obtains 

5L ya (k) = [b 1 + C* r biA(fe)+C p ]5 (1) (fc) + C* 1 ,5< 1) (fe) 
+ [b 1 +Grb 1 A(q 1 ,q 2 ) + C p ] 5 (2) (<h,q 2 ) 

+C v 5f ) {q 1 ,q 2 ) + 



\b2 + \Crb2A{ qi ,q 2 ) 



+ ^Crrb 2 A( qi )A(q 2 ) + ^C pp + C Tp biA{ qi ) 

+b 2 1 CrA(q 1 ) + b 1 C p ] 5 W ( qi )8 W (q 2 ) 
+ [CrvbiA( qi ) + C pv + hC v ] <5 (1) (g 1 )5< 1) (g 2 ) 

(B6) 

- " '^ saidA( qi ,q 2 ) = 



+ lc vv 5l 1) (q 1 )Si 1) (q 2 ), 



where k = q 1 +q 2 , A(q 1 

axctanOgj^iaJA) Where algo 5 Ci)(fe) is rea d as the n = 

111-1-92 I A v ' 

1 term and 8 i - 2 \q 1 ,q 2 ) is read as the n = 2 term of 
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Equation |2o| In the peculiar- velocity effects, 5^ (fc) = 
-ffi 2 rj w (k) and 8 (2 \q 1 ,q 2 ) = -fn 2 2 rj (2) (fc), the terms 
^'(fc) and 7/ 2 '(fc) are the n = 1 and n — 2 terms of 
Equation |21| respectively. In the above expression, fj, is 
the cosine of the angle between the line-of-sight vector, z, 
and the wavevector, fc, i.e., fj, = k-z/k, and / = dl j^"^ is 
the growth rate of structure whose derivative is a function 
of the scale factor, a. 

After we symmetrize the above expression and col- 
lect into first- and second-order expressions, we obtain 
the first- and second-order real-space kernels as 

Z[°\k) = bl (l + Cr ^n(|fc|A) \ +Cp _ f ^ m 



Taking the Fourier transform of both sides, yields the ex- 
act expression relating the real and redshift space quan- 
tities: 

5Lya,s{k) = S Lya (k) (B13) 
+ J d 3 xe-^ (e~ ik * fu * - l) [l + S Lya (x)], 

where u z (k,z) = — ^-r](k,z), and r](k,z) is given by 
Equation[2T] Expanding the exponential as a power series 
and keeping terms up to the second order only, 

<5Ly Q , s (fc) = 6 Lya (k) - 5 v (k) - / d 3 xe- ih - x 



l -blCvvA{ qi )A{q 2 ) + X -C PP + WC P 
+ | [1 + Cr%,? 2 )] - i/C^ + zi?) 
+ \bi (Cr P + 6iCr)[A( 9l )+^( 92 )] 
~^b 1 C v f{n\ + &) - /M12 C v G 2 s) (q x , q 2 ) 



ik z fu z (x)S Lya (x) + ~kzf 2 ul(x) 



(B14) 



Now, inserting Equations [2l] and [31] into |B14| and 
symmetrizing the arguments, we obtain the first- and the 
second-order redshift-space kernels as 



K[ a) (k) = b 1 [l + C T A(k)] + C p + ff/(l-C v ),(B15) 



+ {b 1 [l + CrA(q 1 ,q 2 )]+C p }F^(q 1 ,q 2 ) ^fai.*) = Z^(g u q 2 ) + f^G™ (q u q 2 ) 



A*), 



~^fbiC rv [n 2 2 A(Qt) + £A{q 2 )} 

1 ,2 2 2 r , 

<2^ A'iM2L'i>d- 



(B8) 



+ 2(912/^12/) 



+ 2(912/412/) 



Sfzt\q 2 ) + ^\q 1 ) 
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(B16) 



B2 Redshift-space Lya kernels 

In this subsection, we derive the symmetrized kernels, 
Jfn' , which give the LAE galaxies in redshift space as 



d 3 q 1 /•d 3 q £ _ i 
(2tt) 3 



where 

H12 = (<?! +q 2 ) ■ z/q 12 , 

912 S Iqj+^l, 

<ji z = q x z = q\\X\. 



(B17) 
(B18) 
(B19) 



/n 
d 3 q n 8 D (k-^q i )K i a s) (q 1 ,q 2 ,...,q rl ) APPENDIX C: CALCULATION OF Vb 



x5i(q 1 )6i(q a )...6i(q n ), 



(B9) 



In galaxy redshift surveys we measure positions of 
galaxies in redshift space as opposed to real space. To 
generate the redshift-space expressions for the Lya ker- 
nels we perform the following coordinate transform: 



To calculate the sample variance of the bispectrum (and 
reduced bispectrum), we need the volume of the bispec- 
trum estimator, Vb, used in Equation |48| which deter- 
mines the number of triangular configurations sampled 
at each position. We begin with the definition of Vb'- 



s = x + (l + z)- 



H{z) 



(BIO) 



where s denotes redshift space and v(x) is the line-of- 
sight peculiar velocity. We assume the standard 'plane- 
parallel approximation,' in which the peculiar velocity 
vector and the line-of-sight direction are parallel (or anti- 
parallel), and are chosen to be along the indirection. 
In other words, we ignore curvature of the sky. We can 
rewrite Equation |B11| as 

(Bll) 

■ We can then relate the fluc- 
tuations in real space and redshift space using the mass 
conservation: 



Vb = I d 3 q 1 I d 3 q 2 I d 3 q 3 S D (q 123 j 



Rewriting S D (q 123 ), we find 
d 3 x 



(CI) 



V B = 
Now 



(2,r)3 



i3 ix-q-i 

d q 3 e y3 



qie H1 / d q 2 e H2 I d q 3 e H3 



. , sniffers;) . , .,,,3., 
4?r [ k 3 — Afc + Q(Ak 3 ) ] 



s = x + fu z (x)z, 
where u(x) = (1 + z) 



47T k, 



X 

sin(A;3x) 



Afe, 



(C2) 



[1 + S Lya , s (s)]d 3 s = [1 + S Lya (x)]d 3 x. 



(B12) 



yielding 

Vb = 32 7 rfc 1 fc 2 fc 3 (Afc) 3 f Mkix)sm{k 2 x)sm{k 3 x) ^ 

Jo x 

(03) 
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Evaluating the integral, 

sin(fci:r)sin(fc22:)sin(fc3:r) 



~8 kl 



SJ-i 



sin(fcia;) 



dx 

sm(k 2 x)sm(k 3 x)dx 
lx J sin(A;2a;)sin(fc3a;)da; 



and 

dB 
dln(H) 



dB dki glnfci 
dk~i ainfci d\n{H) 
dB dk 3 <91nfc 3 



+ 



+ 



OB dk 2 dlnk 2 
dk 2 ainfc 2 dln(H) 
dB dm 



dk 3 dlnk 3 dln(H) 
dB dfi 2 . dB 



+ 



d^i dln{H) 



(D2) 



(l/i 



dx 
2^ 



ix(k\i± — k2+ks) 



a ix(k 1 /j,+k 2 - 



_k 2 — k 3 
ki 

k 2 + k 3 
ki 



k 2 — k 3 
fci 

+ 
fci 



<9/i 2 ain(H) 91n(-ff) ' 

where for notational convenience we have dropped the 
dependant variables on all functions. For the determina- 
tion of the actual derivatives, in the evaluation of the 
second-order kernels, we write ki% = £3, fei3 = k 2 , and 
k 23 = fei; and ^12 = —^3, Mi3 = — fJ-2, and /i 2 3 = —Ml- 
Following this we find, for example 



dki 



dki 



dk 3 
dh 



0. 



dki 



dh 



and 



dk 23 
dki 



(D3) 



Here, we have used 
sin(fc2a:)sin(fc3a:) = 



1 ( * 



xk? ixkrt 



+ e 



(C5) 



We apply similar logic when calculating the derivatives 

(C4) 

with respect to k 2 and k 3 as well as when performing the 
derivatives with respect to /ii,/^2, and ^3. 

When performing the derivatives of the bispectrum 
with respect to, e.g., fci, we hold k 2 , k 3 , m-i fi 2 and m 
fixed, despite k 3 and m depending on fci as per the defini- 
tion of the chain rule. We have chosen to write Equations 
|Dl|and|D2|in the form above, to then insert the following 



V B = 87r 2 fc 1 fc 2 fc 3 (Afc) 3 x < 



1 normal triangles 
1 — 6 if fci = fc2 + k 3 
So if k 2 — ki + k 3 or 
k 3 = fci + fc2 



where Q is the Heaviside theta function, and Oo 



dB 



dB dki ainfci dB dk 2 dink 



dln{D, 



+ 



dki Slnfci dln{D A ) 
dB dk 3 ainfc 3 

+ dfc3~dlnfc 3 dln(D A ) 
dB d^ 2 dB 



dk 2 ainfc 2 dln(D A ) 
dB dm 



+ 



dm dln(D A 

dm 



Once we consider the triangular condition, and that 
ki > k 2 > k 3 , we obtain 



derivatives of fci and m from Shoji et al. (20091, 
ainfci 



dln{D A ) 

dlnki 
(C6) dln{H) 

dm 



Mi 



= -Mi 



dln(D A 

and 

dm 
dln(H) 



-fu(i - m ), 



= -M»(! - Mi 



(D4) 



(D5) 



(D6) 



(D7) 



APPENDIX D: DERIVATION OF THE 
BISPECTRUM DISTANCE DERIVATES 

To calculate the Fisher matrix for the angular di- 
ameter distance and the Hubble rate, we need to 
calculate derivatives of the bispectrum with respect 
to the distances, namely, d B ( fc i.fc2^fc3.Mi.M2,M3) anc j 



■fc2.fc3.Mi.M2.M3) 

Sln(D A ) 

dg(fci,fc2.fc;^Mi.M2,M3) _ However, we are dealing with the 
variables fci, k 2 , k 3 , my M2 and ^3, which are not all in- 
dependent of one another. Hence we perform the chain 
rule with respect to the independent functions, fci, 
k 2 , k 3 (6 12 , /ii,<j>), mi M2(mii0)> and ^3 (fei , fe, #12, Mi> <fi)- 
Here, we have written the dependent variables, fc3, m 
and ^3, as functions of the independent variables fci, fc2, 
#12, m and cj>, which fully describe the shape and orien- 
tation of the bispectrum triangles. 
Performing the chain rule, 



APPENDIX E: GENERATION OF THE 
iV-DIMENSIONAL LIKELIHOODS 

In Section [3] we use the the Fisher matrix to compute 
the one and two-dimensional likelihoods marginalizing 
over the remaining parameters in a model of size n. In a 
similar way one can compute iV-dimensional likelihoods 
marginalizing over n — N parameters: 



£?(Xi, x 2 , x N ) 



(El) 



exp 



J2 F ik(Fh 



Fli 



+2^2xi3 

j>i 



ft i 7 



n-JV 

J2 F tk (F k 



-'Fu 



(E2) 



dmdln(D A ) dmdln(D / 



(Dl) 
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